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We investigate the dynamical Casimir effect and its detection with Rydberg atoms. The photons 
are created in a resonant cavity with a plasma mirror of a semiconductor slab which is irradiated 
by periodic laser pulses. The canonical Hamiltonian is derived for the creation and annihilation 
operators showing the explicit time-variation in the couplings, which originates from the external 
configuration such as a nonstationary plasma mirror. The number of created photons is evaluated as 
squeezing from the Heisenberg equations with the Hamiltonian. Then, the detection of the photons 
as the atomic excitations is examined through the atom-field interaction. Some consideration is 
made for a feasible experimental realization with a semiconductor plasma mirror. 



PACS numbers: 42.50.Pq,42.50.Lc,42.50.Ct,32.80.Ee 



I. INTRODUCTION 

The quantum nature of the vacuum provides a vari- 
ety of physically interesting phenomena, including the 
Casimir effect [Ij. The so-called dynamical (nonstation- 
ary) Casimir effect (DCE), as well as the static force, 
has been investigated extensively [3-[23j (and references 
therein), where photons are created from the vacuum 
fluctuation due to nonadiabatic change of the system 
such as vibration of a cavity or expansion of the uni- 
verse. It is, however, difficult experimentally to realize 
the mechanical vibration of the cavity with a sufficient 
magnitude at the resonant frequency ~ GHz which is 
required to create a significant number of photons for 
detection. As a feasible alternative, it has been proposed 
recently that the oscillating wall can be simulated by a 
plasma mirror of a semiconductor slab which is irradiated 
by periodic laser pulses (see also Refs. [H, 03)- 

In this paper, we investigate quantum mechanically 
the photon creation via the DCE and its detection with 
Rydberg atoms. We particularly intend to examine the 
experimental realization of DCE with a plasma mirror of 
a semiconductor slab QjlH!]- In Sec. [TTJ the canonical 
Hamiltonian for DCE is derived in terms of the creation 
and annihilation operators, where the field operators are 
expanded simply with the initial modes. Then, in Sec. 
IIIII the time- varying frequencies and squeezing couplings 
of the Hamiltonian are calculated in an effective 1+1 di- 
mensional scalar field model with a plasma mirror. They 
exhibit the enhancement of effective wall oscillation for 
the DCE which is simulated by the nonstationary plasma 
mirror. In Sec. IIV1 the number of photons created via 
the DCE is evaluated as squeezing from the Heisenberg 
equations for the creation and annihilation operators. 
The results appear to agree essentially with those ob- 
tained by the usual instantaneous-mode approach. In 
Sec. El we investigate the excitation process of Rydberg 
atoms through the atom-field interaction, which is uti- 
lized to detect the created photons. Some conditions on 
the physical parameters are clarified for the efficient pho- 
ton detection. In Sec. IVI1 the experimental realization of 



DCE with a semiconductor plasma mirror is discussed. 
Section IVlII is devoted to a summary. 



II. CANONICAL HAMILTONIAN 

We consider a scalar field in 3+1 space-time dimensions 
as an effective description of the electromagnetic field in 
a resonant cavity. The Lagrangian is given by 

C = \e{^tm 2 ~ ~(V0) 2 - im 2 (x,^ 2 (1) 

(fi = c = 1) 0,1 01, 01, 0. Here, e(x,t) and m 2 (x,i) 
represent the dielectric permittivity and conductivity (ef- 
fective "mass" term), respectively, in the matter region 
such as a semiconductor slab. As specified later, they are 
space-time dependent, simulating the boundary oscilla- 
tion. Conventionally, the instantaneous modes f a (x, t) 
(real, orthonormal and complete) at each time t with 
time- varying frequencies u) a (t) are adopted according to 
the boundary oscillation: 

[-V 2 + m 2 (x,# a (x,() = e (x,^(i)/ a (x,f) (2) 

with the orthonormalization 



e(x,t)/ Q (x,t)/ /3 (x,t)d 3 x = 8 a p/[2uj a {t)]. (3) 



Instead, we here specify the particle representation sim- 
ply in terms of the initial modes 



/°(x) = / o (x,t = 0), LU° a =w Q (i = 0). 



(4) 



The canonical field operators in the Heisenberg picture 
are expanded with the creation and annihilation opera- 
tors a' a (t) and a a (t) as 



(5) 



0(x,t)=5>a(*) + a£ (*)]/(> (x), 
n(x,t)= e (x,0)^^ha a W +ffla (t)]/»(x), (6) 
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where Il(x, t) = d£/d(f) = e(x, t)4>(x, t). Then, the 
canonical Hamiltonian is presented by the usual proce- 
dure as 



H F (t) 



1 r n 2 



+ 0[-V 2 + m 2 (x, } d 3 x 



v 2 U(x,i) 
^0J a (t) (aia a + ]- ) + ^(t)al 



[g af) (t)alal - glp(t)apa a 



(7) 



where the space-integral is taken over the whole region 
V which is fixed suitably (not time-dependent) accord- 
ing to the physical setup, as illustrated later for the case 
of a cavity with a nonstationary plasma mirror. [The 
usual oscillating boundary may also be described as a 
periodic shift of the region of a high potential wall repre- 
sented by m 2 (x, t).} The explicit time-dependence of the 
Hamiltonian Hp(t) in Eq. ([7]) represents the variation 
of the couplings which originates from the nonstationary 
behavior of the c-number external quantities e(x, t) and 
m 2 (x, t). The second-order field equation (Klein-Gordon 
equation) is derived from the Heisenberg equations for 
4>(x, t) and II(x, t). 

The mode frequencies w Q (i), intermode couplings 
H a p{t) and squeezing terms g a p(t) are calculated by con- 
sidering the orthonormality of f® (x) which obey the wave 
equation with e(x, 0) and m 2 (x, 0): 

w a (t) =uj° a + fi aa (t) = u° a + 5uj a (t) , (8) 
^ t p(t)=2G c aP (t)+2G^ p (t), (9) 
9ap(t) = -i[-G% p {t) + G^ p (t)], (10) 

f!^/° (x )/°(x)^, (11) 
sv(t) £A(x,t) 

mi(x,t)/0(x)/°(x)d 3 x. (12) 



SV(t) 



The space-integrals for G^'J'(i) are evaluated actually 
in the subregion SV(t) (C V), possibly time-dependent 
when a moving boundary is considered, where e(x, t) and 
m 2 (x, t) vary in time as 



6 A 1 (x,t) = 6- 1 (x,t)~£- 1 (x,0), 

m^(x,i) = m 2 (x,t) — m 2 (x, 0). 



(13) 
(14) 

Here, G^(0) = with e^foO) = and m^(x,0) = 
at t = 0, as the Hamiltonian Hp(0) is diagonalized in 
terms of the initial modes /„(x). 

Similar formulas are presented for the effective Hamil- 
tonian with the instantaneous modes @, E3. This ef- 
fective Hamiltonian involves even the time-derivatives 
of the mode functions since the quantum time evolu- 
tion is traced along the instantaneous modes. On the 
other hand, in the present approach the time evolution 
is viewed on the initial modes according to the Heisen- 
berg equations. The canonical Hamiltonian is calculated 
without the time-derivatives of the mode functions, and 



applicable readily for various physical setups, e.g., the 
case of a plasma mirror, clarifying its dependence on 
the experimental parameters. There may be some claim 
concerning the ambiguity on the particle representation 
and photon number since the basis modes are changing 
during the DCE. This ambiguity is, however, spurious 
physically (but might be essential for the case of the ex- 
panding universe, which is beyond the present scope). 
In fact, the instantaneous modes return to the initial 
modes at each period of the oscillation, where the pho- 
ton number operators of the respective descriptions co- 
incide with each other by definition. We can check ex- 
plicitly that when the mode functions are not deformed 
largely in time, as usually considered, this canonical 
treatment provides essentially the same result for the 
DCE as the instantaneous-mode approach. The effects 
of the intermode couplings will be less significant in the 
instantaneous-mode approach, where the Hamiltonian is 
diagonalized at each time. Anyway, the intermode cou- 
plings are usually off resonant, providing subleading con- 
tributions to the DCE. 



III. VIBRATION WITH A PLASMA MIRROR 

We next calculate the time-varying frequencies and 
squeezing couplings of the Hamiltonian for DCE in an 
effective 1+1 dimensional scalar field model with a non- 
stationary plasma mirror which is realized with a semi- 
conductor slab irradiated by periodic laser pulses [15]. 

The dielectric response of the plasma is given by 
e(uj) = e%[l — (idp/cu 2 )} with the plasma frequency lu p — 

(n e e 2 /eiin») 1 / 2 in terms of the effective electron mass m* 
and the conduction electron number density n e propor- 
tional to the laser power W^ascr /pulse. This response for 
the dispersion relation, e(w)oj 2 = t\UJ 2 — (n e e 2 /m*), can 
be taken into account in the slab region [I, I + 5] around 
x = I with a thickness <5(<§; L) as 

e(x, t) = ei(t), m 2 (x, t) — m 2 (t) = n e (t)e 2 /m», (15) 

where m 2 (0) = for Wi aser (0) = 0. (The spatial distri- 
bution of the conduction electrons along the x direction 
may also be considered readily.) The instantaneous mode 
functions are given as 



fk(x,t) 



Dsinkx [0,/) 

Be ik'x + c e - ik ' x [I, I + S] : slab 

Asink[x - S + £(t)] (l + 5,L] 



(16) 



with the dispersion relations 



Q 2 = (k 2 + ki)/e = (fc' 2 + k 2 + rn 2 )/ £l (17) 

(fc' = i\k'\ for fc' 2 < with large m 2 ), where is the 
momentum in the orthogonal spatial two dimensions (not 
shown explicitly) [l^, O, Hi}. The Dirichlet boundary 
condition is adopted at x = 0,L with sinfc[i — 5 + £(t)] = 
0, corresponding to the case of TE modes. The case of 
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TM modes can be treated in a similar way by adopting 
m 2 (x,t) = [(dn e /dx)e 2 /(k 2 l m^] 

The diagonal couplings 5u>k(t) and gukif) are specifi- 
cally calculated in Eqs. ([g])-(fI5]) with Eq. ([T6]) for f%(x) 
at t = as 



Suj k (t) 



{t) + 5 m (t)]/L, 
(i/2)u> k [-8 e (t)+S m (t)}/L. 



(18) 
(19) 



Here, the effective wall oscillation is enhanced as 



S e (t)/5 ~ -[ ei (0)/c o ][l - ei(0)/ci(t)] sin 2 kl, (20) 
S m (t)/S ~ ln e (t)e 2 /m*e (LO°) 2 } sin 2 kl. (21) 

This effect is almost proportional to the square 
of the mode function around the slab, [f^(l)] 2 oc 
sin 2 fc;, since j! +S [ft(x)] 2 dx ~ [/£(J)] a <S for k'8 ~ 
[ei(0)/e ] 1/2 ((5/i) < 1 at t = 0. If the slab is placed 
at the boundary x = I = 0, sin 2 kl is replaced with 
(kS) 2 /3 ~ (S/L) 2 <C 1, as observed in Ref. [19| claiming 
that the DCE is suppressed in the TE mode. The signifi- 
cant photon creation, however, can take place even in the 
TE mode if the slab is placed apart from the boundaries 
x = 0, L which are the nodes of f%(x) [H [H|. 

The shift in the instantaneous mode of Eq. (fTo]l 
is determined mainly proportional to 6 to give the fre- 
quency modulation fcfc(t). The diagonal squeezing cou- 
pling <jkk(t) is then calculated with the formulas for the 
effective Hamiltonian (UGH. After some calculations we 
find the relations 

5u) k (t) - Su k (t), g kk (t) ~ [i/2u) k (t)]g k k(t), (22) 

where the change of dielectric is assumed to be small, 
|ei (£) — ei (0) | <C ei(0) as usual [19j. These relations in 
Eq. (|2"2"]l ensure almost the same result for the DCE in 
the canonical and instantaneous-mode approaches (ex- 
cept for the small contribution of the off-resonant inter- 
mode couplings) . This will be checked numerically in the 
next section. 

The above calculations of Sujk(t) and gkk(t) are valid 
up to \SuJk(t)\/cj° k = \S e (t) + S m (t)\/L ~ 0.1, which is 
still a significant enhancement of the effective displace- 
ment |<5 e ,m| 3> S for the DCE. The present approach on 
the fixed basis, however, does not work effectively in an 
extreme situation where the mode functions are largely 
deformed in time with |<5wfc(f)| ~ cj£. In such a case the 
instantaneous-mode approach is rather suitable though 
the deformation of the mode functions cannot be treated 
perturbatively [23]. Anyway, as seen in the following, a 
reasonable deformation to induce \Su>k(t)\/uj1 ~ 0.01—0.1 
is sufficient to create a significant number of photons for 
detection with atoms. 



IV. PHOTON CREATION AS SQUEEZING 

Once the Hamiltonian is presented in terms of the 
creation and annihilation operators, the time evolution 



for the DCE is determined by the Heisenberg equations 
6 a (t) = i[H F (t),a a (t)] and dt(t) = i[H F (t), a+ (t)]. It is 
described as the Bogoliubov transformation, 

a a (t) = A a j3 (t)ap + B* afs (t)oJ , (23) 

al(t)=A* ap (t)al+B afj (t)a p . (24) 

The master equations for the Bogoliubov transformation 
are derived from the Heisenberg equations as 

A a /3 = -iu a (t)A a p - i^ aj (t)A 7/ 3 + 2g ai B 1 p, (25) 
B a p = iuj a (t)B a p + in* aj (t)B 7 p + 2g* ai A lj3 , (26) 

where the intermode couplings are renamed suitably as 

McryCl — da-y) ~ > H-u-y with /J aa = 0. 

In the following, we illustrate the characteristic fea- 
tures of DCE by concentrating on a single resonant 
mode with time- varying frequency uj(t) = loq +5u(t) and 
squeezing coupling g(t) (the mode index u k" omitted). 
The intermode couplings will not provide significant con- 
tributions since they are fairly off resonant generally for 
the nonequidistant frequency differences [H|,[l3|,[2l[. The 
master equations read 

A = -iu(t)A + 2g(t)B, B = iu(t)B + 2g*(t)A (27) 

for the Bogoliubov transformation, 

a{t) = A(t)a + B*(t)a\a)(t) = A*^ + B(t)a.(28) 

The solution is expressed as squeezing and phase rotation 
1, 

A(t) = coshr(t)e t<t>A(t) ,B(t) = sinhr^e^^, (29) 

with the initial condition ^4(0) = 1,5(0) = 0, ensuring 
|A(t)| 2 -|£(i)| 2 = l. 

An analytic solution for A(t) and B(t) is obtained in 
the RWA (rotating-wave approximation) by replacing 

u)(t) — > luo + (8u>) (average), (30) 
g(t) -> (g) ne~ lQt (Fourier component), (31) 

where loq = w(0). By noting the time-evolution of the 
number operator aJ(t)a(t) = \B(t)\ 2 aa^ + . . ., we obtain 
the photon creation via DCE (vacuum squeezing) as 

'|at(t)a(i)|0> = \B(t)\ 2 

smh 2 X t (|A| < \2(g) n \) 



n 7 (£) 



2(.9>s 



X 



\ X \H 2 (|A| = \2{g) n \) (32) 
sin 2 | x |t (|A| > |2( 5 ) n |) 



with the effective squeezing rate 
X = V\2(g^F 



A 2 . 



(33) 



Here, the detunin g A is introduced for the frequency Q 
of laser pulses 12l Il3| as 



Q = 2(lu q + (Sui) + A). (34) 
The resonance condition for DCE is then given by 

fi(resonance) = 2(uj + (Sui)), (35) 
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involving the average shift of the frequency (6u>) [T3 [23| , 
rather than the naive condition O = 2lo . If f2 = 2ujq 
is taken with A = — (Suj), the squeezing rate x is signif- 
icantly reduced, even possibly becomes imaginary with 
n 7 (i) < 1 oscillating as sin 2 \x\t. The photon damping 
with the factor e~ Tt due to the cavity loss should further 
be taken into account, where 

r = oj /Q (36) 

with the cavity quality factor Q. Hence, the threshold 
condition for the squeezing by DCE is placed as 



x > r/2, 



(37) 



which is readily satisfied with a large enough Q. 

We have solved numerically the master equations in 
Eq. (I27p without the KWA. The time-varying couplings 
are taken typically as u>(t) = ujq + (8u) (1 — cosfW) and 
g(t) = 2{g) n (l - cos fit), where |2( ff ) n | - \{Su)\/2 as 
indicated in Eqs. (|T8)) and (jT9]) for the plasma mir- 
ror. The instantaneous- mode solution has also been ob- 
tained by considering the relations Su>{t) = 5ui(t) and 
g(t) = [i/2Cd(t)]g(t) in Eq. fl22J|. In Fig. [TJ the pho- 
ton creation n 7 (t) in the early stage of DCE is plotted 
for iVpuise = i(Q/27r) < 30 (the number of periodic laser 
pulses). The results of the canonical and instantaneous- 
mode approaches are shown with the solid and dotted 
curves, respectively. Here, the parameters are taken typ- 
ically as (5<jj) = 0.02wo, 2(g)n = iO.OlojQ, and A = 
(upper curves), — (Sui) (lower curves) for ft in Eq. (f3~4"|). 
We can see that n 7 (t) increases rapidly via the DCE on 
the resonance with O = 2(ljq + (Soj)) (A = 0), while 
n 7 (t) does not grow for £1 = 2w (A = — {Soj)) due to the 
effective detuning brought by the average shift (6oj) . In 
Fig. [2l the photon creation n 7 (i) is plotted through the 
DCE period for N pulse = t(fl/2n) < 300. The squeezing 
rate is determined from this plot to be x — O.Olwoj a $ 
indicated in Eq. (133)) with A = 0. This result confirms 
that a large number of photons can be created via the 
DCE with a reasonable squeezing rate x ~ O.Olwo when 
the laser pulses are applied many times. It is also found 
that the canonical and instantaneous-mode approaches 
provide almost the same result (except for the small con- 
tribution of the off- resonant intermode couplings). The 
analytic solution under the RWA in Eq. (|3"2"|) overlaps 
almost with the instantaneous-mode result though it is 
not plotted explicitly in Figs. [T]and[2j 

We briefly discuss the effect of the intermode couplings. 
Specifically, the coupling ii\2a\a,2 + ^i\ 2 a \ a i between the 
modes 1 and 2 becomes resonant under a condition = 
3w? -> UJ2 ~ w i = w ft for the case of the TE m 
and TE115 modes in a cubic cavity due to the relation 
(l 2 + l 2 + 5 2 ) 1 / 2 = 3(1 2 + l 2 + l 2 ) 1 / 2 . Then, through 
this resonant intermode coupling the significant photon 
creation occurs in both the modes 1 and 2 as n~,i(t) ~ 
n 7 2 (t) , increasing the total of photon numbers [13|, [2l], 
|23I |. The photons of the mode 2 are, however, fairly off 
resonant with the Rydberg atoms tuned to detect the 
photons of the mode 1. Hence, they cannot be detected 
efficiently. 
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FIG. 1. Photon creation n 7 (t) (linear plot) in the early stage 
of DCE for iVpuise = t(Q/2n) < 30 (the number of periodic 
laser pulses). The results of the canonical and instantaneous- 
mode approaches are shown with the solid and dotted curves, 
respectively. The parameters are taken typically as {Su) = 
0.02^0, 2(g)n = iO.Olwo, and A = (on-resonance: upper 
curves), —{Slu) (off-resonance: lower curves) for Q. 
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FIG. 2. Photon creation n 7 (t) (log plot) through the DCE pe- 
riod for iVpuise = t(Q/2n) < 300. The results of the canonical 
and instantaneous-mode approaches are shown with the solid 
and dotted curves, respectively. The parameters are taken 
typically as (Slo) — 0.02o;o, 2{g)n = zO.OIcjo, and Q. — 2.04uo 
(A = 0). 



V. DETECTION WITH RYDBERG ATOMS 

The photons created via the DCE are detected suit- 
ably by Rydberg atoms with principal quantum number 
n fa 100 and transition frequency ~ GHz jlH l23l|. Ryd- 
berg atoms may be treated as a two-level system with a 
transition frequency uj e for the resonant photon absorp- 
tion with oj e ps ujq. They are initially prepared in the 
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lower level \g), and injected into the cavity. A part of 
these atoms are excited to the upper level |e) by absorb- 
ing the photons, and detected outside the cavity as the 
signal of photons. Recently, a high-sensitivity measure- 
ment of blackbody radiation has been performed at a fre- 
quency 2.527 GHz and low temperatures 67 mK - 1 K by 
employing a Rydberg-atom cavity detector with a newly 
developed selective field ionization scheme for n w 100 
(the atoms excited by absorbing photons are selectively 
ionized by applying an electric field) [2~4j . Here, we note 
that in order to observe purely the vacuum squeezing via 
DCE, the cavity should be cooled well below 100 mK 
to suppress the thermal photons as n 7 (thermal) ^ 1. 
In fact, if photons are present initially with an expecta- 
tion value (a^a), they are also amplified by the DCE as 
(l + 2\B(t)\^a). ' 

Consider that N Ryd Rydberg atoms (actually N Ryd ~ 
100 — 1000 [3]), which are all prepared at the lower 
level \g), are injected into the cavity to detect the cre- 
ated photons after the period of DCE, for simplicity of 
argument. (The following features for the photon de- 
tection are essentially valid even if the atomic beam is 
injected continuously during and after the DCE, as dis- 
cussed later.) The n 7 photons and N Ryd atoms (all lo- 
cated at the same position for simplicity) are coupled 
with the Jaynes-Cummings Hamiltonian under the RWA 
as 

H AF = n^N R J d ~{aD + + ^D-). (38) 

(The effect of the counter-rotating terms is negligible 
near the resonance.) Here, the collective atomic spin-like 
operators are defined (in the Schrodinger picture) 25] by 

D + = £ \e)(g\ (i) /y/N^, (39) 

i=l 
N Kyd 

D_=J^\g){e\ (i) /y/N^, (40) 

i=l 

and the complex phase for k is absorbed in the atomic 
levels. The single atom-photon coupling n is explicitly 
given by 

K = d v 4V2^y(|/ (x 1 )|/|/°(x )|) (41) 

in terms of the magnitude of the electric dipole transition 
matrix element d, the cavity volume V and the mode 
function /°(x), where Xi and Xo represent the atomic 
position and the antinode, respectively. The collective 
atom-photon coupling is suitably defined by 

R = n^N Ryd . (42) 

The single atom-field coupling is typically k~3x 10 3 s~ 1 
at the antinode for the Rydberg atom of principal quan- 
tum number n w 100 with ui e « uip ~ 1.5 x 10 10 s _1 
(2.4GHz x 2tt) and V ~ (0.1m) 3 [H,|25j]. Then, the col- 
lective coupling amounts to R ~ 10 5 s _1 ~ 10 _5 Wo f° r 
N Ryd ~ 10 3 , which is still much smaller than the reso- 
nant frequency uj e k, loq. 



The commutation relations among the collective oper- 
ators are given by 

[D+,D-] = D Z 

N Ryd 

= E[l e >< e l(*)-|3>(slM]/^Ryd, (43) 

i=l 

[D z ,D±] = ±(2/N Ryd )D±. (44) 

The operators N e and N g to represent the populations of 
the upper and lower levels \e) and \g), respectively, are 
given by 

Ne=Yl l e )( e lw = (N Ryd /2)(l + D z ), (45) 
i=i 

JVRyd 

N g = £ \g)(9\(i) = (N Ryd /2)(l - D z ), (46) 

i=l 

satisfying the completeness 

JVRyd 

N e + N g =Y, Oe)(e| (i) + \g)(g\ (l) ] = N Ryd . (47) 

i=l 

The created photons are detected by counting the num- 
ber of excited atoms which is represented by N e with 
eigenvalues 0,1, . . . , N Ryd . The initial atomic state is 
prepared as 

|0 e ) = \gi,g 2 ,---,gN Ryd ), (48) 

which is an eigenstate of N e with zero atomic excitation 
satisfying D_|0 e ) = 0. The one-excitation state is gener- 
ated as 

|le) = £>+|0 e > 

and so on for the multi-excitation states. 

The Heisenberg equations are derived by taking 
the total Hamiltonian Ha + Haf + Hp with Ha = 



{N Ryd /2)uj e D z for the free atomic system: 

a = —iojQa — iRD - , (50) 

D_ = —iuj e D- + iRaD z , (51) 

D z = -i(2/N Ryd )R(aD+ - a)D_). (52) 



We solve these equations perturbatively to see the evolu- 
tion of the atomic excitation N e {t) — (N e {t)). First, Eqs. 
(ISTR) and ([51]) for a(t) and D_(t) = D\_(t) are integrated 
up to the first order of R with the initial atomic operators 
D±(ti) in Eqs. (|3"9")l and (|4"0)) and the photon operator 
a(ti) at t = t\ after the DCE with one sequence of N pn \ BC 
laser pulses for the duration 

h = N pulsc {2ir/n). (53) 
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Then, the results are applied to Eq. (l52j) to obtain D z (t) 
up to the second order of R. This determines the atomic 
excitation as 

N e (t) = (N e (t)) = (N Ryd /2)[l + (D z (t))} 

~n 7 (2S/A e ) 2 sin 2 [A e (t-t 1 )/2], (54) 

where the atomic detuning is given by 

A e = u> e — ujq. (55) 

In these calculations, the following relations are con- 
sidered: {a,al}D x + {£>+,£>_} = 2(a t aD z + Z>+D_), 
[a,a^}D z - [£>+,£>_] = 0, (0 e |I>±(ti)|0e) = 0, 
<0 B |£»+(ii)D_(ti)|0 e ) = 0, (0 e |U,(ti)|0 e ) = -1, and 
(0|a t (£i)a(£i)|0) = n 7 (the photons created via the 
DCE). Note here that N e (t) < N Ryd with (D z ) w -1 in 
the early epoch of photon detection (the linear regime). 
Although it is difficult in practice to trace exactly the 
time evolution beyond the linear regime for the system 
of the many atoms interacting with the resonant cavity 
mode, we may survey the essential features for the atomic 
excitation to detect the photons as follows. 

Suppose that n 7 ^> N Ryd , namely the photons are cre- 
ated much more than the Rydberg atoms, as desired and 
feasible experimentally. Then, the atomic excitation is 
eventually saturated as N e (t) <~ («;i) 2 n 7 ~ Nji yd for 
t l/(/Cy/n^), which is expected by extrapolating Eq. 

(|54|) roughly up to Rt ~ ^/AR y d/n 7 <C 1 near the reso- 
nance A c « (henceforth t — t\ — > t). This excitation 
process may be viewed as the onset of Rabi oscillation 
between \g) and |e) at a rate 

which takes place almost independently for the N Ryd 
atoms in the presence of the large field (many photons 
with n 7 > N Ryd ). 

On the other hand, if n 7 < N Ryd though less inter- 
esting experimentally, the excitation is exchanged be- 
tween the atoms and field as N e (t) ~ n 7 /2 on average 
for Rt ~ 1. This may be understood from the fact that 
the interaction Hamiltonian i?AF in Eq. (|38j) describes 
the oscillation with a rate Q e ~ k = fsy A^Ryd between 
the atomic and field operators in the linear regime. The 
collective atomic excitation can be treated as a quantum 
oscillator, satisfying approximately the bosonic commu- 
tation relation [D-,D+] ~ —(D z ) « 1 with n 7 <C AR y d 
in Eq. (|43p . that is D + and D_ act as the creation and 
annihilation operators, respectively (25| . 

The cavity loss eventually becomes significant for t > 
1/r. Then, the atomic excitation is also relaxed with a 
rate 

f 4(^/r) 2 r (s < r/4) 
ie \ r/2 («>r/4) 1 j 

through the transition |e) — > Ig) + 7 and the loss of the 
emitted photon in the cavity [25J. We also note that the 



atom-field interaction terminates when the atoms transit 
through the cavity. The atomic transit time is given by 

t tT = L/v = I^ 1 , (58) 

where v and L are the atomic velocity and the cavity 
length, respectively. We have typically 

r tr ~^=3xl0V\ (59) 
0.1m 

which is comparable to the single atom-field coupling k. 
By considering these damping effects, we realize that the 
created photons are detected efficiently with the atoms 
under the conditions, 

r> e >r,r tr , (60) 

Ttr > T e . (61) 

That is, the atomic excitation should take place for t ~ 
fi" 1 before the significant loss of the created photons due 
to the cavity damping (r > 2r e ), and the actual cutoff of 
the atom-field interaction by the atomic transit (r tr ). It 
is also required that the excitation damping (r e ) induced 
by the cavity loss does not become significant before the 
atoms transit through the cavity (r tr ). 

As investigated so far, if the photons are created copi- 
ously via the DCE with n 1 3> AT Ryd , they are detected 
by the atomic excitation as 

N e (t tl ) ~ N Ryd /2. (62) 

Here, the condition Q e > T tr is less restrictive, requiring 
merely n 7 > (r tr /K) 2 ~ 1 for r tr ~ K. The atomic 
detuning may be suppressed readily as A e < fl e , e.g., for 
fi e -3x lO^" 1 with k~3x lO^ 1 and n 7 - 10 6 . The 
conditions fi e > T and r tr > T e ~ 2{R/T) 2 T (R < r/4) 
imply lower and upper bounds, respectively, on the cavity 
quality factor, 

(cV0/v^ £ Q Z (u Q / K ){r tr / K )/N Ryd , (63) 

where ujo/k ~ 5 x 10 6 . These bounds are combined as a 
requirement for the number of created photons, 

n- ( > ( K /r tr ) 2 iV 2 yd » iV Ryd . (64) 

For example, we estimate Q ~ 5 x 10 3 and n 7 ~ 10 6 
for r tr ~ K and N Ryd ~ 10 3 . This range of Q meets 
consistently the condition R < r/4 for T e . 

On the other hand, if T e = T/2 (R > r/4) the condi- 
tion Ttr > r e places a significant bound 

Q>m/r tr ~5x 10 6 . (65) 

This range of Q meets consistently the condition R > r/4 
for r e . We also note that N e (tt r ) ~ n 7 /2 for n 7 < N Ryd - 
In this case with £l e ~ R, the condition f2 e > T implies 
R > r/4. Hence, the above range of Q in Eq. (j6"B"|) is 
effective either for n 7 > N Ryd or n 1 < N Ryd . 

The atomic beam may be injected continuously 
through the period of DCE. Then, we can show that 
the atomic excitation is squeezed together as N e (t) ~ 



7 



{n/(jjf)) 2 n 1 {t) during the DCE. This atomic excitation 
is usually smaller than ^V Ryd ~ 100 — 1000, e.g., for 



k/uj 



10 5 and n-v < 10 



10 



Anyway, the created 



photons are detected with the atoms efficiently after the 
DCE. 



VI. EXPERIMENTAL REALIZATION 

We now discuss a feasible experimental realization of 
DCE with a semiconductor plasma mirror [l5|,[23[. Based 
on the analyses presented so far for the DCE and pho- 
ton detection, we can find desired values for the physical 
parameters. 

The photons are created as 



1 

r 



2xtl 



i e 27r( x /a;o)iV puls 



(66) 



with the squeezing rate x f° r the resonant mode, where 
t\ = -^ P uise(27r/f2) and f2 ~ 2wo (see also Fig. [2|. Hence, 
the desired number n 7 of created photons places a re- 
quirement for the squeezing rate as 



as 

X = \2(g) a \ ~ 0.01wo(r n /0.1)(Ww/0.01/iJ). (70) 

Here, the factor m represents the Fourier component 
(g)ae~' int of g(t), which may be optimized by suitably 
designing the time-profile Wiaser(t) of laser pulse. As 
seen in Eqs. and (|55|) . the tuning of fi is required 

for the resonance by taking into account the average shift 
{6u))/w ~ 0.01 -0.1. 

As for the photon detection, the analyses in Sec. |V] 
indicate that roughly Nnyd/2 ~ 100 atomic excitations 
are detected per mean atomic transit time itr ~ 0.1ms 
for the creation of ri 7 ~ 10 6 — 10 s photons via the DCE. 
The quality factor of cavity should be chosen suitably 
to ensure the efficient atomic excitation and detection. 
Specifically, Q ~ 5 x 10 3 in Eq. JB3J or Q > 5 x 10 6 
in Eq. (|6"51) . We note that even if an excessive amount 
of photons (n 7 ^> 10 s ) are created, their detection is 
actually limited by the number of Rydberg atoms N^y^ ~ 
100—1000. After the detection, the photons remaining in 
the cavity are relaxed finally as n 1 — >• for t > 10ms 3> 
T ,tt T ', namely the field returns to the vacuum. Then, 
the subsequent rounds of photon creation and detection 
are performed repeatedly. 



ln(4n 7 ) 

27r7Vpulsc 



(67) 



Typically, x ~ 0.01w to obtain n 7 ~ 10 6 - 10 s with 
-^Vpuise = 300 laser pulses, where the threshold condition 
X > r/2 for the DCE is also satisfied sufficiently with 

Q £ io 3 - 

The effective displacement in Eq. (|21[) is achieved by 
applying a laser power Wiaser/pulse for the period T = 
2ir/n ~ 0.2ns: 



5 m /L ~ (n s e 2 /e m >t )L/7r 2 



(68) 



where sin 2 kl = 1 for dcfinitcncss (the slab is placed in 
the middle of cavity I = L/2), luqL ~ tt, and n s = n e 5 
(ex Wiaser) is the surface number density of electrons. We 
may readily obtain (n s e 2 /eom*)L ~ 1 with a reasonable 
laser power Wi ascr /pulse ~ O.Ol^iJ/pulse [2J], achieving 
a significant displacement S m ~ 0.1L. In this case, the 
conductivity effect S m in Eq. (f2Tj) dominates over the 
dielectric effect S e in Eq. for ei(0) — 1 - 10 and 

£i(0) < | ei (t) j [the photon damping by the complex e%(t) 
does not exceed the squeezing by the DCE mainly with 
<5m] ■ We estimate the variation of the mode frequency as 

5u ~ (S m /L)u> ~ 0.1wo(Wi aser /0.OlMJ). (69) 

By noting the relation \Su\ ~ \2g\, the desired squeezing 
rate for the DCE can be obtained in Eq. f|33[) with A = 



VII. SUMMARY 

We have investigated quantum mechanically the pho- 
ton creation via DCE and its detection with Rydberg 
atoms, specifically considering the experimental realiza- 
tion in a resonant cavity with a plasma mirror of a semi- 
conductor slab irradiated by laser pulses. The canoni- 
cal Hamiltonian for the DCE is derived in terms of the 
creation and annihilation operators showing the explicit 
time- variation which originates from the external config- 
uration such as the nonstationary plasma mirror. Then, 
the photon creation is evaluated as squeezing from the 
Heisenberg equations. This confirms that a sufficiently 
large number of photons can be created via the DCE 
with a reasonable squeezing rate when the laser pulses 
are applied many times. The atomic excitation process 
to detect the photons is described with the atom-field in- 
teraction, which clarifies the conditions for the efficient 
detection. Based on these analyses, desired values of the 
physical parameters are considered for a feasible exper- 
iment for DCE and its detection with a plasma mirror 
and Rydberg atoms. 
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